ON THE Z 2 -COHOMOLOGY OF 3D POLYHEDRAL 
APPROXIMATIONS 



ROCIO GONZALEZ-DIAZ, JAVIER LAMAR, AND RONALD UMBLE 



Abstract. The cohomology algebra of an object encodes important geometri- 
cal information not captured by its cohomology groups. Let / = (Z 3 , 26, 6, B) 
be a 3D digital image, let Q(I) be the associated cubical complex and let 
dQ(I) be the subcomplex of Q(I) whose maximal cells are the quadrangles 
of Q(I) shared by a voxel of B (the foreground) and a voxel of 1? x B (the 
background). We show how to simplify the combinatorial structure of dQ(I) 
and obtain a polyhedral complex P(I), which is homeomorphic to Q(I) but 
with fewer cells. We provide an algorithm that computes cup products on 
H* (P(I);1,2) directly from the combinatorics. The computational method in- 
troduced here can be effectively applied to any polyhedral approximation of a 
3D object. 



1. Introduction 

This paper completes the work proposed in our prequcl-abstract [B] by provid- 
ing additional insights, examples, proofs and a corrected version of formula ([1} in 
Theorem 14. II All modules herein are assumed to have Z2 coefficients. 

Let X be a cellular complex embedded in 3-dimensional space and constructed 
by gluing 3-dimensional polyhedra together along common faces (see [S])- At the 
vector space level, the cellular cohomology H*(X) is generated by the connected 
components of X, homotopy classes of non-contractible loops modulo commutator, 
i.e., tti(X)/ [ni(X), tti(X)], and surfaces bounding the cavities in X. At the alge- 
bra level, the cup product encodes certain relationships among these generators, and 
improves our ability to distinguish between objects. For example, H* (S 1 V S 1 V S 2 ) 
and H*(S 1 x S 1 ) are isomorphic as vector spaces but not as algebras since 2- 
dimensional cup products vanish in the wedge but not in the product; consequently, 
S 1 V S 1 V S 2 and S 1 x S 1 have quite different topological properties. 

To date, the cup product has seen limited application to problems in 3D image 
processing. In [7j, Gonzalez-Diaz and Real used their 14-adjacency algorithm and 
the standard formulation in [T5] to compute cup products on the simplicial complex 
K (I) associated with a given digital image /. In contrast, Kaczynski, Mischaikow 
and Mrozek [3] showed how to compute the homology of a digital image directly 
from the pixels thought of as cubical complex. Their algorithm, which applies tech- 
niques from linear algebra, demonstrates that the homology groups of a cubical 
space are computable. More recently, Gonzalez-Diaz, Jimenez and Mcdrano [5] 
introduced a method for computing cup products directly from a cubical approxi- 
mation Q(I) of a given digital image / (no additional subdivisions are necessary). 
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Motivated by problems in high-dimensional data analysis, Kaczynski and Mrozek 
[lOj give an algorithm for computing the cohomology ring of a cubical complex of 
arbitrary dimension. The implementation of the method presented in their paper is 
currently in progress. For a geometrical interpretation of cohomology in the context 
of digital images, we refer the reader to [4]. 

A problem that frequently arises in 3D image processing is to efficiently encode 
the boundary surface of a given digital object as a set of voxels. The most popular 
approach to this problem uses a triangulation. While triangles are combinatorially 
simple, and visualization of triangulated surfaces is supported by existing hardware 
and software, the number of triangles required is often large and the computational 
analysis is correspondingly slow. It is desirable, therefore, to consider combinato- 
rial approximations with simpler combinatorics and a subsequent improvement in 
computational efficiency. 

Given a triangulation, we introduce a technique for merging adjacent triangles 
to form polygonal faces of a combinatorially simpler polyhedral approximation. 
Using polyhedra to approximate 3D objects is well-studied problem in the field of 
Computational Geometry (see [UEJG], for example). 

In [llj , Kravatz constructed a combinatorial diagonal approximation on general 
polygons and used it to compute cup products on closed compact orientable surfaces 
thought of as identification spaces of 2?i-gons. However, applications of Kravatz's 
diagonal in more general settings are limited by its dependence on a particular 
indexing of the vertices. Here we introduce a more general and computationally 
effective diagonal approximation, which does not depend upon the indexing. 

The computational method introduced here can be effectively applied to any 
polyhedral approximation of a 3D object. Our method maximizes computational 
efficiency by minimizing the number of cells in a given polyhedral approximation. 
A Matlab implementation of the method is available on the Internet for testing. 

The paper is organized as follows: A review of the standard definitions from 
Algebraic Topology appears in Section [2] and the notions of cup product and di- 
agonal approximation are reviewed in Section [3] Our main result, which appears 
in Section [4] as Theorem l4.1[ gives an explicit formula for computing cup products 
on the polygons in a given polyhedral approximation P of an object embedded 
in 3-space. This is sufficient for computing the cohomology algebra H*(P) since 
H 3 (P) = 0. Let / = (Z 3 , 26, 6, B) be a 3D digital image, let Q(I) be the associated 
cubical complex and let dQ(I) be the subcomplex of Q(I) whose maximal cells arc 
all the quadrangles of Q(I) shared by a voxel of B (the foreground) and a voxel of 
Z 3 \ B (the background). In Section [5] we show how to simplify the combinatorial 
structure of dQ(I) to obtain a polyhedral complex P(I), which is homcomorphic 
to Q(I) but with fewer cells. In Section [6] we produce an algorithm for computing 
cup products on H*(P(I)) directly from the given combinatorics and analyze its 
computational complexity. The algorithm presented here is adapted from the al- 
gorithm for computing cup products in [7|. Conclusions and some ideas for future 
work are presented in Section [7] 

2. Standard definitions from Algebraic Topology 

Wc begin with a review of some standard definitions from Algebraic Topology 
(for details see [13]). 
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Let S = {Sq} q be a graded set. A q-chain of S is an element of S q ; finite 
formal sums of g-chains define additive abelian group structures on S q and S. 
The group of g-chains is denoted by C q (S); the chain group of S is the graded 
group C*(S) = {C q (S)} q . A chain complex of S consists of a chain group C*(S) 
together with a square zero homomorphism d = {d q : C q (S) — > C q -i(S)} q , called 
the boundary operator. For example, consider a triangle Vj, Vk) with vertices 
Vi < Vj < Vk' The boundary of the triangle is the formal sum of its edges, that is, 
d 2 (vi,Vj,Vk) = (vi,Vj) + (vj,v k ) + (vi,v k ). 

Intuitively, a cellular decomposition of a 3D space X embedded in R 3 is a rep- 
resentation of A as a finite union of vertices (0-cells), edges (l-cells), polygons 
(2-cells), and polyhedra (3-cells), which have been glued together in such a way 
that the non-empty intersection of two cells is a cell. A fc-cell is also referred to as 
a k-face. In this paper, the term cellular complex refers to a 3D space X embedded 
in R 3 together with a cellular decomposition. For the precise definition of a cellular 
complex, which is more subtle than one might expect, see [8]. An i-ccll a' G X 
is a facet of a cell a £ X if a is an (i + l)-cell and a' is a face of a. A maximal 
cell of X is not a facet of any cell of A. A cubical complex Q is a cellular complex 
whose 2-cells are squares (or quadrangles) and whose 3-cells are cubes. Note that 
if a cube is in Q, its bounding quadrangles are in Q; if a quadrangle is in Q, its 
bounding edges are in Q; and if an edge is in Q, its endpoints are in Q. 

The g-dimensional cells of a cell complex X provide an important example of 
a graded set. The cellular chains of X, denoted by C* (A) , is the chain group 
generated by the cells of X. The cellular chain complex of X is the chain complex 
(C* (A) , d) , where d is the linear extension of the cellular boundary. A g-chain 
o G Cq(S) is a q-cycle if d q (a) = 0; it is a q-boundary if there is a (g + l)-chain 
fi such that d q +\((i) = a. Referring to the triangle (vi,Vj,Vk) above, the 1-chain 
(7 = (vi,Vj) + (vj,Vk) + (vi,Vk) is both a 1-cycle and a 1-boundary since d\{a) = 
and a = di (vi , Vj , Ufe) . Two g-cycles a and a' are homologous if a+a' is a g-boundary. 
Let Z q {S) and B q {S) denote the g-cycles and g-boundaries of S, respectively. Then 
B q (S) C Zq(S) since a 2 = 0. The quotient H q (S) = Z q (S)/B q (S) is the q th 
homology group of S and the graded group H*(S) = {H q (S)} q is the homology of 
S. An element of H q (S) is a class [a] := a + B q (S). The cycle a is a representative 
cycle of the class [a]. The rank of H q (S) is called the i t?l Betti number and is 
denoted by bi. 

Let (C*(5),9) and (C*(5'),d') be chain complexes. A homomorphism / = {/ g : 
C q (S) —> Cq(S')} q is a chain map if f q -id q = d' q f q for all q. Let f,g : C* (S) —> 
C* (5")} be chain maps. A chain homotopy from f to g is a homomorphism = 
{</> g : Cq(S) -» C 9+ i(S")} g such that 4> q -xd q + d' q+l 4> q = f q + g q for all g. A chain 
contraction of (C*(S), 9) io (C*(S"), 9') is a tuple (/, g, 0, (S, <9), (5", 5')) consisting 
of chain maps / : C*(5) — >• C7*(S") and 5 : C*(S") — > C*(5) such that = Id ,( S n, 
and a chain homotopy : C*(S) — > C*(S r ) from Id c< , (S) to g/. 

Definition 2.1. Let A be a cell complex and let e be a p-cell of A contained in 
the boundary of exactly two (p + l)-cells E and B. Let A' be the cell complex 
obtained from X by replacing E and B with the single cell B' = (E U B) \ int(e), 
and define 

a W — I d(E + B) , x = B'. 
Then and -B have merged into £>' along e. 
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0-cochain 


!>i} 


1-cochain 


{ei,e 4 } 


l-coboundary 


H v i} = {e3,e4} 


1-cocycle 


c = {ei,e 2 } 


1-cocycle 


d = {ei,e2,e3, e^} 


homologous cocycles 


c and <i; since d = c + (5{ui} 



Figure 1 . Example of cochain, cocycle and coboundary. 

Proposition 2.2. If cell complexes X and X' are related as in Definition \2.U then 
H*(X) w H,(X'). 

Proof. There is a chain contraction (f,g,4>, (C*(X),d), (C*(X'),d')) given by 
. /(e) = e + 8(E), f{B) = B', f(E) = 0, f(x) = x for x ± e, E, B; 

• g(B') = E + B, g{x) = x for x ^ B'; and 

• 0(e) = E, 4>(x) = for x ^ e. 

Since is a homotopy equivalence, the conclusion follows. □ 

Cochain groups are linear duals of chain groups. Given a chain complex (C*(S), 
d), a q-cochain is an element of C q (S) := Hom(C q (S), Z2). If we index the q-cells in 
a cellular complex S from 1 to n 9 , their corresponding duals generate C q (S). Thus 
a g-cochain is a Z2-linear combination of these duals, and can be thought of as a bit 
string of length n q . The cochain group of S is the graded group C*(S) = {C q (S)} q . 
The coboundary operator is the homomorphism S = {S q : C q (S) — > C q+1 (S)} q 
defined on a q-cochain c by 5 q (c) = cd q +i; then S 2 = and the cochain complex 
associated with (C*(S),d) is the pair (C*(S),6). 

A q-cochain c is a q- cocycle if 6 q (c) — 0. A q-cochain b is a q- coboundary if 
there exists a (q — l)-cochain c such that 6 = (5 9_1 (c). Two q-cocycles c and c' 
are cohomologous if c + c' is a q-coboundary. The q-cocycles and q-coboundarics 
of S arc denoted by and Bi(S), respectively. Then B q (S) C Z«(5) since 

5 2 = 0. The quotient ^9(5) = Z q (S)/B q {S) is the q" 1 cohomology group of S and 
the graded group H*(S) = {H q (S)} q is the cohomology of 5. An element of H q (S) 
is a class [a] := a + B q (S). 

Remark 2.3. Since the ground ring Z2 is a field, the homology and cohomology 
of a graded set S are isomorphic and torsion free (see [8]). Furthermore, given a 
chain contraction (f,g,(j),(S,d),(S',d')), there is a map <9__/ : C* (S") -> C* (5) 
that sends a chain cr to the cochain d a f defined by 

a tf \ / 1; if a is a non-zero summand of / (//) ; 
OvJ W ~ j 0, otherwise. 

Thus if cr 7^ 0, the cochain c^/ is supported on the subspace generated by all 
/i G C* (5) such that cr is a summand of / (/x) . 

An AT-model for a chain complex (C*(5), 9) consists of a chain complex (C*(if), 
9' = 0) together with a chain contraction (f,g,cj),(C*(S),d),(C*(H),d')). An 
AT-model for (C*(5),c5) always exists and can be computed in C(m 3 ), where 
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m = #5 (see [7] for details). A particular AT- model for (C*(S),d) is denoted 
by (f,g,cf>,(S,d),H). 

Proposition 2.4. Let (/, g, <f>, (S, d), H) be an AT -model for (C*(S),d). 

i. If o G H qi then g q (cr) represents a class in H q (C* {S)). 

ii. The map H q — > H q (S) given by a t— s> [9(c)] extends to an isomorphism 
C q (H)^H q (S). 

Hi. The map H q — > H q (S) given by a t— s> [d a f] extends to an isomorphism 
d-f : C q (H)^Hi (S). 

Proposition 2.5. Let (f,g,4>,(S,d),(S',d')) be a chain contraction. 

i- If(fi,9i,4>i, (S,d),H) is an AT-model for (C*(S),d), then (fig, fgi, ffag, 
(S',d'),H) is an AT-model for (C*(S"), <9')- 

ii. If (f 2 ,g 2 ,(f> 2 ,(S',d'),H') is an AT-model for (C*(S'),d'), then (f 2 f,gg 2 , 
+ 5< /) 2 /, (S, 9), H') is an AT-model for (C*(S), d). 

Proof, (i) Observe that fid(n) = and dgi(a) = for all /1 G C*(S) and a G 
C*(if). Hence figfgi = fi{(f>d + dcf> + U c ,(s))gi = Ziffi = Mc,(ff) and f4>igd + 
dffag = f(<j>id + d(j>i)g = /(Idc.(s) +51/1)5 = ^c(s') +fgifi9- 

(ii) Similarly, observe that f 2 d(n) = and dg 2 (a) = for all /i G C*(S") and a G 
C*(H'). Hence f 2 fgg2 = /2(Idc,(S'))92 = /2#2 = Idc.(ff') and ((/) + g<f> 2 f)d+d(<f> + 
gfof) = (j)d+d<j)+g(4> 2 d+d(j)2)f = <f>d+d<f>+g(Id c ,(s>) +g 2 f 2 )f = Id c „(S) +352/2/- 

□ 

3. Cup Products and Diagonal Approximations 

In this section we review the definitions and related notions we need to compute 
the cup product on the cellular cohomology H* (A) of a cell complex A. 

Definition 3.1. The geometric diagonal on X is the map A : X — > X x X 
defined by A(x) = (x,x). A diagonal approximation on A is a cellular map 
A x : A — > A x Xsuch that 

i. A x is homotopic to A. 

ii. if e is a cell of A, then A x (e) C e x e. 

Hi. A x extends to a chain map C*(A) — s- C*(X) <E) C*(A), in which case the 
geometric boundary d extends to a coderivation of A x , i.e., 

A x d = (Id®<9 + <9® Id) A*. 

A diagonal approximation A x always exists by the Cellular Approximation Theo- 
rem (see [8]). If A x is coassociative, (C*(X), d, A x ) is a differential graded coalge- 
bra (DGC) . 

Example 3.2. The Alexander- Whitney (A-W) diagonal approximation A s on the 
n-simplex s n — (0,1, ... ,n), n > 0, is defined by 

n 

A a <0,l J ...,n> = ^(0,...i)®(i,...n) 

i=0 

and extends to a strictly coassociative, homotopy cocommutative chain map A s : 
C*(s„) — > C7*(s„) <g) C*(s„). Thus (C* (s„) ,d,A s ) is a DGC for each n (see [8]). 
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The 0-cube is the single point 0, the 1-cube is the unit interval I = [0,1], and 
the n-cube is the n-fold Cartesian product I™ = I x • • • x I. A cell of I™ is an n-fold 
Cartesian product u% ■ ■ • u n £ {0, 1, 1} x ™. 

Example 3.3. The Serre diagonal approximation Ai on the n-cube I" is given by 

Ai (0) = <8> and Ai(F) = (0 <g> I + 1 <g> l) xn for n > 1, 
where a®bxc®d = ac® bd. In particular, 

Ai (I 2 ) =(0®I + I®1)x(0®I + I®1)=00®II + 0I®I1+I0®1I + II®11. 
Let 0' = I and I' = 1; then for n > 1 we obtain 

Ai(I™) = ^2 u\ ■••u n ®u' l ■■■u' n . 

u 1 ---u n G{0,l} xn 

Again, Ai extends to a strictly coassociative, homotopy cocommutative chain map 
Ai : C*(P) -> C*(I") ® C*(I n ) and (C, (I n ),d,Aj) is a DGC for each n (see [15]). 

Let A x : C*(X) — > C*(X) Cg> C*(A) be a diagonal approximation and let 
{f,g,<j), (A, d), if) be an AT-model for (C*(X),d). Then C* (if) = if*(X) sw H*(X) 
by Remark 12.31 and Proposition 12.41 Given classes a £ if p (X) and b £ if 9 (X) , 
there exist unique chains a = (d-f)^ 1 (a) £ C p (if) and /3 = (d-f)^ 1 (6) £ C g (if) 
such that a — [d a f] and b — [dpf]. The cup product of representative cocycles <9 Q / 
and dpf is the (p + g)-cocycle defined on 7 £ C p + g (X) by 

(9 a / fy/) (7) = • (d a f ® ^/)A X ( 7 ), 

where " • " denotes multiplication in Z2, and the cttp product of classes a = [<9 a /] £ 
HP(X) and 6 = [dpf] £ if«(X) is the class 

a^ x b= [d a f - x dpf] £ HP+i{X). 

The cup product on if* (A) is a bilinear operation with unit. If A x is homotopy 
coassociative, (if* (X) , -— - x ) is a graded algebra with unity; if A x is homotopy 
cocommutative, (if* (X) , -~ ' X ) is (graded) commutative. In particular, ^— - s and 
^ 7 are commutative products with units. 

Definition 3.4. Let A and X' be cell complexes, let A x : C*(A) C*(A)®C*(A) 
be a diagonal approximation, let (/, g, (j>, (A, 9), (A', 9')) be a chain contraction, and 
let (/',#', 0', (X',d'),H') be an AT-modcl for C*(A'). The induced coproduct 
on A' is the chain map 

A*' := (/ O /)(A x ) ff : C* (A') -> C.(A') ® C.(X'). 

Proposition 3.5. Under the conditions of Definition \3J^ the induced coproduct 
A x i is homotopy cocommutative whenever A x is homotopy cocommutative. 

Proof. Let t : A® A — > A<E>A denote the twisting isomorphism given by r(a<8>6) = 
b ® a, and choose a chain homotopy h : C* (A) — > C* (A) ® C* (A) such that 

r A x + A x = (d <g> Id + Id <g)3) /i + ft<9. 
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Then (/ <g> /) hg is a chain homotopy from tA x > to A X ' as the following calculation 
demonstrates: 

tA x , + A x ,=t (/ <g> /) A x g + (/ ® /) A x5 

= (/®/) (rA x +A x ) ff 

= (/®/)((0®Id + Id®fl)/i + ft3)s 

= (9 ® Id + Id ®a) (/ (g) /) hg + (/ <8> /) fcflia. 

□ 

In practice, we compute "cup products" in homology via an AT-model (/, g, <j>, 
(X,d),H) for C*(X). Having obtained a basis for C*{H) H*(X), we define the 
homology cup product as follows: 

Definition 3.6. Let X be a cell complex, let A x : C*(X) -> C*(X) ® C*(X) be a 
diagonal approximation, and let (/, g, (f>, (X, d),H) an AT-model for C* (X). Given 
a £ ij p and (3 £ define 

a := ]T {• (0„/ ® ^/)(A X ). 9 ( 7 )} 7 £ C*(# P+9 )- 

-yeH p+g 

Proposition 3.7. Let X and X' be cell complexes, let A x : C*(X) — > C*(X) (g> 
C*(X) be a diagonal approximation, let (f,g,(j>, (X,d), (X',d')) be a chain contrac- 
tion, and let (f',g',<j>', {X',d'),H') be an AT-model for C*(X'). Then for a £ H' p 
and (3 £ H' q we have 

a x / f3 = a p. 

Proof. Since (/'/, gg', fi+gftf, (X, d), H') is an AT-model for C7* (X) by Proposition 
[231 part (ii), and C*(H*{X) « H*(X') by Remark [23] we have 

a E {• (5a/'®^/')(A x ,)ff'(7)}7 

= E {' (^/'®^/')(/®/)(Ax) 55 '(7)}7 

= E {■ (d a f'f®dpff)(A x )gg'(>y)} 1 = a~ x l3. 

□ 

4. Cup Products on 3D Polyhedral Approximations 

Traditionally, one uses the standard formulas in [13[ 115] to compute cup prod- 
ucts on a simplicial or cubical complex. In this section, we show how to compute 
cup products on a 3-dimcnsional polyhedral approximation P of a cell complex 
embedded in M 3 . 

Arbitrarily label the n vertices of P with the integers {l,...,n} and represent 
a fc-gon p C P by an ordered fc-tuple of integers p = (ii, . . . ,ik), where i\ = 
min{ii, . . . ,ik}, i\ is adjacent to ik, and ij is adjacent to ij+i for 1 < j < k. The 
statement in Theorem 14.11 will use the following notation: The symbol i m nA := 
max {%2, ■ ■ ■ , ik}] the Z2-coefficient Xj = if and only if ij < ij+i; and edges of 
interest are denoted by { Uj = {h,ij)} 2 < j < k and {ej = (i J ,ij+i}} 2 < j < k -i ■ 



8 



ROCIO GONZALEZ-DIAZ, JAVIER LAMAR, AND RONALD UMBLE 



Theorem 4.1. Let P be a polyhedral approximation. The A-W diagonal induces a 
coproduct A P on each polygon p = (ii, . . . , £&) C P given by 



m(k)-l 

A P (p) =(ii) ®p+p® (i OT (fe)> + ^ (w 2 + e 2 H h e 3 -_i + Aj-ej) ® ej 

fe-1 

(1) + [0- + *j)ej+ej+i + --- + e k -i+u k ]®e j . 

j=m(k) 

Proof. Consider the triangulation {Tj_i = ij, ij+i)} 2 <j<k-i °^ P anc ^ n °t e that 
the A-W diagonal on Tj-\ is given by 



A s (Tj_ x ) = Aj ((ii) (8 Tj_i + Tj-x ® (ij) + u j+1 ® e,) 
(2) + (1 + Aj) ((ii) <g> Tj_i + Tj-i ® + Uj ® e,-) . 



Our strategy is to inductively merge these triangles until p is recovered and the 
induced coproduct is obtained. We proceed by induction on j. 

When j = 2, set p\ = T\ and note that either i 2 < 13, in which case J m (3) = *3 
and A 2 = 0, or i 2 > i 3 , in which case i m ($\ = i 2 and A 2 = 1. If i m (v\ = 13, formula 
(UJ) gives the non-primitive terms [u- 2 + A 2 e 2 ] ® e 2 . Since A 2 = this expression 
reduces to m 2 ® e 2 = A 2 U3 <8> 63 + (1 + A 2 ) u 2 ® e 2 . On the other hand, if i m ($) = *2, 
formula (|T|) gives the non-primitive terms [(1 + A 2 ) e 2 + 143] ® 63. Since A 2 = 1 this 
expression reduces to M3 <§) e 3 = A 2 it3 €5 e3 + (1 + A 2 ) u 2 ® e 2 . In either case, formula 
(JT|) agrees with the A-W diagonal on p\ . 

Now assume that for some j > 3, formula ([T]) holds on pj„ 2 = (ii, , ij) . Merge 

Pj_ 2 and Tj_i along and obtain Pj-i = (ii, . . . , ij+i) ; we claim that formula ([T]) 
also holds on Pj-i- Recall that 



A P (pj-i) = (/»/) A s5 = (/ g /) A s fo_ 2 + T^x) . 



Either ij + i > i m u), in which case Xj = and i m (j+i) = b'+ii or l 'm(j) > m 
which case £ m (j+i) = im(j)- First assume that ij+i > i m (j)- Following the proof of 
Proposition 3.3, define / (uj) = u 2 + e 2 + • • • + e^-i, / (pj-2) = 0, and / (Tj_i) = 
Then 



ON THE Z2-COHOMOLOGY OF 3D POLYHEDRAL APPROXIMATIONS 



Apfe-i) = (/® ®Tj_i +Tj_i ® (ij+i) ®ej- 

m(j)-l 

+ (ii)®Pj-2+JJj-2 8>(*m(j))+ X! ( u 2 + e 2 H h A s e s ) 8 e s 

s=2 

j-1 

+ [(1 + A s ) e a + e a +i + • • • + ej-_i + itj] (8) e s } 

s=-m(j) 

= ® (ij+i) 

m(j)-l 

+ (u 2 + e 2 + • • • + ej-i) <g ej + (u 2 + e 2 + • • • + A s e s ) <g> e s 

s=2 

+ V] [(1 + A s ) e s + e s+ i H + ej-i + (U2 + e 2 + h e^-i)] <8> 1 

s=m(j) 

j 

= (h) ®p 3 -i +Pj-i ® (ij+i) + X] (" 2 + e 2 H 1- A s e s ) g> e s , 

s=2 

which verifies formula ([T]) in this case. On the other hand, if i m (j) > define 
/ = e i + M j+i> / (Pj-2) = Pj-ii and / ( T j-i) = 0- Then 

Apfe-i) = (/®/){A i ((i 1 )<g)Tj_i + r i _ 1 (8)(^)+u j+1 (8)ej) 
+ (1 + Aj) ((n) <g> + Tj-i (8 (ij+i) + ® e 3 -) 

to(j)-1 

+ (ii) <8>Pj-2 +Pj-2 ® («m(j)) + ( u 2 + e 2 H hA s e s )®e s 

s=2 

J-1 

+ ^2 [(1 + A s ) e s + e s+ i H h ej_i + (8 e s } 

s=m(j) 

= (n) a Pj-i +P3-1 ® (*m(j)) + AjMj+i ® ej 

m(j')-l 

+ (1 + Aj) (ej + Mj+i) ® + (u 2 + e 2 H h A s e s ) ® e s 

s=2 

j'-i 

+ ^2 [(l + ^) e s + e s+H h e 3 -_i + (ej + itj+i)] ® e s 

s=m(j) 

m(j)-l 

= (ii)<8>Pj-i+Pj-i®(i m (j))+ (w 2 + e 2 H hA s e s )®e s 

s=2 

i 

+ ^2 [(1 + A s ) e s + e s+ i H h e 3 - + Uj+i] <g> e s , 

s=m(j) 

which verifies formula (JlJ in this case as well and completes the proof. □ 
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Figure 2. Left: A digital image I = (Z 3 , 26, 6, B); the set B 
consists of 8 unit cubes (voxels). Right: The quadrangles of dQ(I). 




Figure 3. Critical configurations (i), (ii) and (hi) (modulo reflec- 
tions and rotations). 

5. 3D Digital Pictures and Cellular Complexes 

Let 7 be a 3D digital image and let Q(I) be an associated cubical complex. 
In this section we introduce a simplification procedure, which produces a cellular 
complex P{I) homeomorphic to dQ(I) with significantly fewer cells. 

Consider a 3D binary digital picture / = (Z 3 , 26, 6, B), where Z 3 is the underlying 
grid and B (the foreground) is a finite set of points of the grid fixing the 26- 
adjacency for the points of B and the 6-adjacency for the points of Z 3 \ B (the 
background). The cells of Q{I) are unit cubes centered at the points of B with 
faces parallel to the coordinate planes (called the voxels of I), together with their 
quadrangles, edges, and vertices. The subcomplcx dQ(I) consists of all cells of Q(I) 
that are facets of exactly one (maximal) cell of Q{I), and their faces. Note that 
the maximal cells of dQ(I) are all the quadrangles of Q(I) shared by a voxel of B 
and a voxel of Z 3 \ B (see Figure [5]) . 

Wc perform a simplification process in dQ(I) to produce a cellular complex P(I) 
homeomorphic to dQ(I) such that the maximal cells of P{I) arc polygons. But first, 
we need a definition. 

Definition 5.1. A vertex v G dQ(I) is critical if any of the following conditions 
is satisfied: 

i. v is a face of some edge e shared by four cubes, exactly two of which intersect 
along e and lie in Q(I) (see cubes w\ and in Figure [3]). 

ii. v is shared by eight cubes, exactly two of which are are corner-adjacent and 
contained in Q{I) (see cubes si and S2 in Figure [3]). 
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(a) lb) 
Figure 4. (a) N v <- 92, 93,^4}; (b) facets oip<- {ei, e 2 , e 3 , e 4 , e 5 , e 6 , e 7 , e 8 }. 




xz_337_sch4.tif xz_338_sch4.tif xz_339_sch4.tif xz_340_sch4.tif xz_341_sch4.tif xz_342_sch4.tif xz_343_sch4.tif 




_344_sch4.tif xz_345_sch4.tif xz_346_sch4.tif xz_347_sch4.tif xz_348_sch4.tif xz_349_sch4.tif 

Figure 5. A jiMRl of a trabecular bone. 

in. v is shared by eight cubes, exactly two of which are corner-adjacent and 
not contained in Q(I) (cubes t\ and ti in Figure [3]). 



A non-critical vertex of dQ[T) lies in a neighborhood of dQ(I) homeomorphic to 
R 2 (see [12]). 

Algorithm 15.21 presented below processes the non-critical vertices of dQ(I) to 
obtain the cellular complex P{I). Initially, P(I) = dQ{I). For a vertex v £ P{I), 
let N v be the set of 2-cclls q £ P(I) incident to the vertex v. If N v defines a region 
R v homeomorphic to a disc, then N v is replaced by a new 2-cell p in P(I), which 
is the union of the cells of N v . The edges of P(I) incident to v and the vertex 
v are removed from P(I) (see Figure SJ. Observe that the maximal cells of the 
final cellular complex P(I) are polygons and P(I) has fewer cells than dQ(I). Wc 
need some terminating conditions. For example: (1) terminate when the number 
of edges of the polygons in P(I) reach some specified maximum; or (2) terminate 
after merging the set N v of coplanar 2-cells of P(I). Example 15.31 demonstrates the 
differences that arise from these different terminating conditions. 

Algorithm 5.2. Algorithm to obtain the cellular complex P{I). 

Input: the cubical complex Q(I). 
Initially, P(I) := dQ(I); 

V := ordered set of non-critical vertices of dQ(I). 
While 3v £ V such that R v is homeomorphic to a disc do 

remove v from P(I) and V ; 

remove the edges incident to v from P(I) ; 

remove the 2-cells of N v from P(I); 
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Figure 7. (a) The cellular complex P(^I) for 10 edges as upper 
bound on p 6 P(/j,I). (b) The cellular complex preserving 
geometry. 



add a new 2-cell p to P(I) 

which is the union of the 2-cells of N v . 
Output: The cellular complex P(I) . 

Observe that Algorithm 15.21 uses the ordering on the set of non-critical vertices 
V C dQ(I) to select the next non-critical vertex. To the best of our knowledge, this 
is the first algorithm to appear that produces a cellular complex with polygonal 
maximal cells by removing non-critical vertices. 

Example 5.3. Let /i/ be a /J.MRI of a trabecular bone of size: 85 x 85 x 13 voxels 
(see Figure [5] in which /il is given by a sequence of 13 2D digital images of size 
85 x 85). The number of quadrangles in dQ(fil) is 23136 (see Figure After 
applying Algorithm 15.21 to dQ(fil) with the terminating condition 1 (the number 
of the edges of the polygons in P(fil) is at most 10), the number of polygons of 
P(fJ,I) is 1044 (see Figure[7]a). If we only consider the set of coplanar polygons N v , 
then the number of polygons of P((J>I) after applying Algorithm 15.21 is 11895 (see 
Figure 0b). 
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30 



Figure 8. Representative 1-cycles on the cellular complex P(^il). 

Now, consider the output P(I) of Algorithm 15.21 The chain complex 
(C*(P(J)), d) can be described as follows: 

• The chain group C q (P(I)) is generated by the g-cells of P(I). 

• The value of d q : C q (P(I)) — > C q -\{P(I)) on a g-cell is the sum of its facets. 

• The boundary of a sum of g-cells is the sum of their boundaries. 

Example 5.4. Consider the cellular complex P(^I) displayed in Figurc[7Jb, which 
was produced by an application of Algorithm 15.21 to a /iMRI of a trabecular bone 
of size 85 x 85 x 10 voxels. Table [1] exhibits the homology obtained by computing 
an AT-modcl for P(fil). These results count the number of connected components, 
holes and cavities. Representative 1-cycles are shown in Figure |U 

Table 1. Results of the homology groups computation for the 
cellular complex P(/J,I) shown in Figure [TJb. (see Figure [8]). 



Cellular Complex b bi b 2 
P(jjJ) 10 9 17 



6. Computing the Cohomology Algebra H*(P(Ij) 

The following algorithm computes cup products in H*(P(I)) in terms of A P . 
Since H 3 (P(I)) = 0, all non-trivial cup products in H*(P(I)) arc products of 
1-cocycles. 



14 



ROCIO GONZALEZ-DIAZ, JAVIER LAMAR, AND RONALD UMBLE 




Figure 9. Representative 1-cocycles on the cellular complex P{(xl). 




Figure 10. Left: cell complex A. Right: cell complex B. 



Algorithm 6.1. Computing Cup Products on P(I). 

INPUT: A cubical complex Q(I) associated to a 3D digital image /. 
Compute a polyhedral cell complex P(I) homeomorphic to dQ{I) 

using Algorithm 15 . 21 . 
Compute an AT -model (f,g,4>,(P(I),d),H) for (C*(P(I)),d) 

using the algorithm given in H7J . 
Initially let A be a bl ( b i +1 ) x hi zero matrix. 

Let Hi = {ai,..., a bl } and H 2 = {jx, . . . , 7& 2 } . 
For i = 1 to i>i do 

For j = i to b\ do 

For k = 1 to 62 do 

A(i + j, k) = -(d a J ® d aj f)(A P )g( lk ); 
Output: the matrix A. 
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Observe that 

b 2 

}]{A(i + j, k)}j k = a t P otj 
fc=i 

by Definition 13.61 and the matrix A is symmetric since A P is homotopy cocommu- 
tative by Proposition [3751 

Let to be the number of voxels of /. Since the number of cells of Q(I) is at most 
27 ■ to, i.e., 0(m), the number of critical vertices of dQ(I) is 0(m). The overall 
computational complexity of Algorithm 6.1 is 0(m + ?i 3 + 62 • n ■ k ■ b\), where bi is 
the i th Betti number; its three components are determined as follows: 

• The computational complexity of Algorithm 15 . 2 1 to obtain the cellular com- 
plex P(I) is 0(m) since, in the worst case, the number of edges in P(I) 
incident to a non-critical vertex v S P(I) is 6 and the number of 2-cell of 
N v is 12. 

• Let n be the number of cells of P(I) and k the number of vertices of the 
biggest polygon of P(I). The computational complexity of the algorithm 
given in [7J to obtain an AT-model (f,g,(j>,(P(I),d),H) for (C*(P(/)),9) 
is 0{n 3 ). 

• The complexity to compute a row of A is at most 0(n ■ k 2 ■ b\), since for a 
fixed 7 S H2, 3(7) has at most n summands, A p has k summands and for 
a 1-cell a £ P{I), and f(a) has at most b\ summands. 

Finally, observe that, in most cases, k « n « to and bi « n « m for £ = 1,2. 
Then, the overall complexity in most cases is 0(n 3 ). 

Example 6.2. Starting from the results obtained in Example 15.41 and applying 
the formula given in Theorem 14. 11 representative 1-cocycles are shown in Figure [9] 
Table [2] shows the cup product. 



Table 2. Results of the computation of the cup product on the 
cellular complex shown in Figure [71b. Only non-null results are 
shown. 





(c*2 , CI4) 


(«3 , a s ) 


(«5 , 07) 


(«6 , a 9 ) 


(«7 , ag) 


716 


1 


1 


1 


1 


1 



Table [3] illustrates the dramatic improvement in computational efficiency realized 
by removing faces and non-critical vertices. 

Example 6.3. Let DBrain be a 3D digital image obtained from a MR imaging 
of the brain. After applying Algorithm 15.21 with terminating condition (2), i.e., 
terminate after merging the set N v of coplanar 2-cells of dQ(DBrain), we obtain 
a cellular complex P (DBrain) with 2433 polygons (see Figure [TT]Left). After 
computing an AT-model (/, g, <j>, (P(DBrain)), d, H) for (C*(P(DBrain)),d) we 
obtain 60 = 2, b\ = 19 and 62 = 2. Representative 1-cycles are shown in Figure 
[TTjRight. Table |4] shows the resulting cup product. 
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Table 3. Time (in seconds) and the results of the computation of 
the cup product on the cellular complexes shown in Figure 1101 



Cell complex Number of 2-cclls Time to compute the cup product 

A 1638 28.00 sec. 

B 46 1.04 sec. 

(ai,a 2 ) (oi,o 3 ) (01,04) (o 2 ,o 3 ) (a 2 ,o 4 ) (03,04) 

7 1 1 1 



Table 4. Results of the computation of the cup product for the 





cellular complex 
shown. 


shown 


in Figure 


fm Only non- 


null results 


are 


72 


1 


o 2 o 3 

1 


O2Q5 
1 


Q2O17 
1 


a 3 a 5 
1 


03 o 6 
1 


0:307 
1 


o 3 o 9 
1 


72 


0:30:13 
1 


0506 
1 


05«7 
1 


0509 
1 


050:13 
1 


06012 
1 


6 0i3 

1 


06016 
1 


72 


OsOig 
1 


09017 
1 


O11O12 
1 


012019 
1 


013015 
1 


01401! 
1 


j Oi6Ql9 
1 


018019 
1 
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Figure 11. Left: cell complex P(DBrain). Right: Representa- 
tive 1-cycles. 



7. Conclusions and Plans for Future Work 

Given a 3D digital image /, we have formulated the cup product on the coho- 
mology of the cellular complex P{I) obtained by simplifying the cubical complex 
dQ(I). The algorithm proposed here is valid for any encoding of a 3D digital object 
given as a set of polyhedra. Our ultimate goal is to compute cup products on a 
general n-dimcnsional cellular complex over a general ring directly from the given 
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combinatorial structure (without subdivisions). Our strategy is to apply some stan- 
dard geometric constructions such as forming quotients, taking Cartesian products, 
and merging cells. 
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